Effect of impurities on quasi-2D quantum antiferromagnet 
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We have studied the static and dynamic properties of quasi-two-dimensional quantum antiferro- 
magnets (AF) diluted with spinless impurities using spin-wave theory and T-matrix approximation. 
We show that the spectrum of a 2D AF at long wavelengths is overdamped at arbitrary concentra- 
tion of spinless impurities. The scattering leads to a new length scale i/a ^ e'^^*^, x being impurity 
concentration and a the lattice spacing, beyond which the influence of impurities on the spectrum 
is dominant. Although the dynamical properties are significantly modified we show that 2D is not 
the lower critical dimension for this problem. Thus, in low-dimensional systems with disorder the 
connection between static and dynamic quantities is not straightforward. Our results are in quan- 
titative agreement with the recent Monte Carlo simulations and experimental data for S = 1/2, 
5=1, and S = 5/2. We have also proposed experiments which can further test the results of our 
theory. 
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I. INTRODUCTION 

The discovery of superconducting cuprates has at- 
tracted much attention to the properties of diluted 2D, 
QHAF [p]-p^, the subject intensively studied some 30 
years ago in the context of magnetism in diluted mag- 
netic alloys ||ll|;|l2| . Traditional view of the effect of local 
disorder on the spectrum of an ordered 3D antiferromag- 
net is that at long wavelengths the form of the spectrum 
is not modified. The only effects are the reduction of 
hydrodynamic parameters and a weak damping. 

In this work we study the problem of impurities in 
2D QHAF within the linear spin-wave theory using the 
T-matrix approach. The spin-wave Green's function is 
evaluated by summing all multiple-scattering diagrams 
that involve single impurity. We find that the scatter- 
ing leads to a new length scale i/a ^ e^/^^ beyond 
which the influence of impurities on the spectrum is dom- 
inant. We associate this length scale with the localization 
length of spin excitations. We show that the dynamical 
structure factor 5{k, w) for a^^ » A; » consists of 
three parts (we use units h — ks = (j)a broadened 
quasiparticle peak with a width given by Tk — xc^k, 
Co = 2^/2SJa; (ii) a non-Lorentian localization peak at 
u! = loq cq£~^, (Hi) a flat background of states between 
Lu = cok and uj = luq. Thus for every k-state some weight 
is spread from the high energies to the low energies. For 
k ^ the quasiparticle and localization peaks merge 
into a broad incoherent peak that disperses in momentum 
space. We calculate the static magnetic properties and 
flnd a quantitative agreement with both MC simulations 
and experimental data. We show that at T = the stag- 
gered magnetization is given by M{x, 0) « 5* — A — Bx, 
the factor A « 0.2 is the contribution of the zero-point 
fluctuations of the spins and B ~ 0.21. We find that 
Tn{x)/Tn{0) ~ 1- As X where As = n-2/TT+B/{S-A). 
This result gives A1/2 — 3.2 and — 2.6 which agree 



very well with experiments up to x ~ 0.25. 

The systems discussed in this paper are modeled by 
the site-diluted quantum Hciscnbcrg antiferromagnet: 



(1) 



where pi = 1 (0) if i's site is occupied (unoccupied) by 
the spin S. We focus on the problem of tetragonal or 
square lattices with in-plane, J, and out-of-plane, Jj_, 
nearest-neighbor exchange constants, (ij) denotes sum- 
mation over bonds. In the systems of interest J ^ J±. 
We begin with the Hamiltonian (|l]) in the linear spin- 
wave approximation which is split into the pure host and 
impurity part and, after the Bogolyubov transformation, 
is given by (in the units of 45 J) 
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with two-component vectors A 
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(2) 
(3) 

and 



2x2 scattering potential matrices Vk,k' (for details see 
]l^), I runs over the impurity sites. We are interested 
in the Green's function of the Hamiltonian (^) modi- 
fied by random impurity potentials (^. The Green's 
function is a 2 x 2 matrix defined in a standard way 
and the T-matrix equation for the Hamiltonian (0) is 



given by: 



-V, 



l.j-L 

k,k' 



T.y^'Gl{u)T'^>i,{u:), 



with / = -^{B), partial waves are restricted to in-plane 
/i — s,pa,d harmonics and G'q(w) is the bare Green's 
function. The T-matrix equations can be solved and the 
subsequent averaging over random distribution of impuri- 
ties transforms T-matrix into the spin-wave self-energies 



Sk(w) = a:(5k-k' I], 
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Summation 
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of the Dyson-Belyaev diagrammatic series for the Green's 
functions gives G = & + G^SG and then the transverse 
component of the neutron scattering dynamical struc- 
ture factor is proportional to the linear combination of 
the spectral functions 

5^(k,c.)=/(k)5H[AiiM+ Al^{cj) + 2Al'iu;)] , 

where A'^{lu) — — ilmG^ (a;), kinematic form-factor 
/(k) = ttS* (l-7k)/wk, gioj) = [l + UBiu})], and n_B(cj) = 



1] ^ is the Bose distribution fmiction. 



The static properties of the system are calculated from 
the spin- wave expression of the averaged on-site magnetic 
moment: 



\{s:)\ 



where bosonic averages can be expressed through the 
spectral functions. 



vicinity of a "localization peak", w « wq, where the spec- 
tral function rises sharply from the shallow background 
states '--^ X to a peak of the height ^ 1/x and then van- 
ishes in a singular fashion as uj approaches zero. The 
spectral function Al^{uj) is shown in Fig. ^ for a repre- 
sentative value of impurity concentration x = 0.1. One 
can clearly see the features we have discussed above: the 
broadened quasiparticle peak, the localization peak, and 
the states between them. As the k decreases all the men- 
tioned structures merge. 

The average on-site magnetic moment for randomly di- 
luted AF with the averaging over magnetic sites M{x) — 
J2i \Si\/Nm can be expressed through the integral of the 
spectral functions as: 



Mix,T) = S 
6M{x, T) 
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SM{x,T) , 



(4) 



II. RESULTS 



where A — X^k'^k — 0-1966 is the zero-point spin devia- 
tion, subscript R denotes retarded. 
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FIG. 1. The spectral function A]^{uj) for the wave- vectors 
k = 0.005, 0.01, and 0.02 along the (1, 1) direction, k = 0.005 
is of the order of uq. 
the picture. 



(ui) for each k is normalized to fit 



A detailed analysis of the spectral function gives the fol- 
lowing picture. At the wave-vectors much larger than 
^0 (f-^k 3> cjo), that is at the wavelengths shorter than 
a characteristic length £ ^ e^'^/*'^, the spectral function 
has three distinct regions in uj. First, is a vicinity of a 
quasiparticle peak, uj w Wk, where the spectrum has a 
regular Lorentzian form with the pole at (Dk and width 
7k given by the perturbative result. Second, the interme- 
diate range of energies, loq < uj 'ti oj-k, where the spectral 
function is independent of uj and corresponds to an al- 
most flat, shallow (~ x) background of states. Third, the 
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FIG. 2. Average staggered magnetization v.s. x. Our re- 
sults from Eq. (^) (solid line) , Monte Carlo data (circles, Ref. 

NQR data (diamonds, Ref. |l||), and the fit of Monte 
Carlo data from Ref. (dotted line) are shown. 

Numerical integration of the expression in Eq. (Q) 
gives the suppression rate of the staggered magnetization 
M{x) ~ M(0) - Bx with B = 0.209(8). For 5* = 1/2 it 
gives the slope of the normalized staggered magnetiza- 
tion M{x)/M{Q) ~ 1 - Bx/{S - A) ~ 1 - 0.691(5) • x. 
The staggered magnetization v.s. x is presented in Fig. 
I for S* = 1/2. The recent Monte Carlo data Refs. §,0, 
and NQR data Ref. are also shown. One can see a 
very good agreement of our results with numerical data 
up to high concentrations. 
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FIG. 3. Tjv(a;)/rjv(0) v.s. x for S = 1/2. Our results (solid 
line), /iSR (diamonds) magnetic susceptibility (circles) Q, 
and ESR (squares) M data. 



any anomaly associated with the spectrum and are fi- 
nite up to the concentration close to the classical per- 
colation threshold. These results are in a quantitative 
agreement with the NQR, /iSR, ESR, and magnetic sus- 
ceptibility measurements in different compounds as well 
as with the Monte Carlo data. Thus, in low-dimensional 
systems with disorder the connection between static and 
dynamic quantities is not straightforward. Altogether 
this provides a self-consistent picture of the effects of dis- 
order in low-dimensional quantum antiferromagnets. 
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At T > the staggered magnetic moment consists of 
the quantum, T = 0, and thermal, T-dependent, parts. 
For the true 2D system at x = and T > thermal 
fluctuation destroy the long-range order which manifests 
itself as a log-divergency of the thermal correction to the 
magnetization. The 3D coupling provides a cut-off to this 
divergency in a quasi-2D problem which yields the finite 
value of the thermal correction and the finite value of the 
Neel temperature whose mean-field value can be found 
from the condition M(x, Tn) = = S - A- 6M{x, T/v). 
Suppression rate of the Neel temperature can be readily 

obtained ~ 1 - As x 1 ~ x(^Tr - ^ + . For 

S = 1/2 this gives Ai/2 = 3.196(5) and for 5* = 5/2 it is 
A5/2 = 2.600(4). We plot our results for Tn{x)/Tn{Q) for 
the case of 5' = 1/2 in Fig. || together with experimental 
data. One can see that our results agree very well with 
the experiments up to a rather high doping level x ~ 0.25. 



III. CONCLUSIONS 

We have studied the problem of diluted 2D and quasi- 
2D quantum Heisenberg antiferromagnets in a tetrago- 
nal lattice making use of linear spin-wave theory and 
T-matrix approach. We have shown that the spin-wave 
spectrum is strongly modified by disorder which indicates 
magnon localization on a length scale i, exponentially 
large in 1/x. This new length-scale appears explicitly in 
the dynamic properties such as the dynamical structure 
factor 5(k, uj) which can be measured directly in neutron 
scattering experiments, and the magnon density of states 
N{uj), which is directly related to the magnetic specific 
heat. The measurement of such quantities will provide 
a direct test of our theory. Furthermore, we show that 
the static properties such as the zero-temperature stag- 
gered magnetization and Neel temperature do not show 
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